Tight Binding Molecular Dynamics Studies of Boron Assisted Nanotube Growth 



E 

Institut de Ciencia de Materials de Barcelona - 



Hernandez * and P. Ordejon 

CSIC, Campus de la Universitat Autdnoma de Barcelona, 08193 Bellaterra, 
Barcelona, Spain 



o 
o 
o 

(N 

G 
G 



Bergische Universitat - Gesamthochschule Wuppertal, 



I. Boustani 

FB 9 - Theoretische Chemie, 
Germany 



Gauflstrafle 20, D-42097 Wuppertal, 



A. Rubio and J. A. Alonso 

Departamento de Fisica Tedrica, Universidad de Valladolid, 47011 Valladolid, Spain 

(February 1, 2008) 

In this paper we report a theoretical study of the effects of the presence of boron in growing car- 
bon nanotubes. We employ a well established Tight Binding model to describe the interactions 
responsible for the energetics of these systems, combined with the Molecular Dynamics simulation 
technique and Structural Relaxation calculations. We find, in agreement with the previous theoret- 
ical/experimental work of Blase et al. [Phys. Rev. Lett. 83, 5078 (1999)], that boron favors (n,0) 
(zig-zag) tubular structures over (n,n) (arm-chair) ones by stabilizing the zig-zag edge. Further- 
more, it is shown that boron has the effect of delaying the tube closure process, a fact which could 
explain the improved aspect ratio experimentally observed in nanotubes synthesized in the presence 
of boron. Our dynamical simulations lead us to propose a mechanism through which this extension 
of the closure time can be explained. 



G 

O 

o 



> 
o 

CO 
CN 

\o 
o 
o 
o 

-I— > 



i 

G 

O 

o 



X 



I. INTRODUCTION 

The discovery of carbon nanotubes by Iijima in R 
has marked the starting point of a scientific revolution 
This discovery has opened a whole new perspective in the 
nanoscopic regime. aLMaterials Science and Engineering. 
Their mechanicalHJ, electricat3'E£l and magnet idllTO 
properties provide ample opportunity for the fabrication 
of nano-scale devices. Indeed some such-jde vices have al- 
ready been reported in the literatureE2rE3. Since Iijima's, 
discovery nanotubes of other chemical composition have 
also beerij^aynthesized, such as B x C y N z composite nan- 
otubesEj - E3, the so-called inorganic nanotubes consisting 
of layers of M0S2 or WS2CJ Ej, or MO2 nanotubesE3. 
All these compounds have phases which consist of lay- 
ered structures, and this has lead to the prediction that 
other materials also capable of crystallizing in layered 
structures can in principle produce nanotubes. Indeed, 
theoretical arguments have been presented in_the liteca- 
tiKe for tlw.viabilijte of BNEU, BC 3 El, BC2NB, GaNO, 
BE3, GaSecJ and FO nanotubes. Interestingly, nanotubu- 
lar structures can also be constructed from biochemical 
compounts, such as peptides, as demonstrated experi- 
mentally by Ghadiri et aLcil and theoretically by Carloni 
et alH. 

The tubes first detected by Iijimatl were multi-wall 
nanotubes (MWCNT's), i.e. concentric shells of cylin- 
drical shape, in which each shell is separated from the 
next by approximately the same distance as the inter- 
layer spacing in graphite. Each shell can be character- 
ized by a pair of indices, (n,m), which determine how the 
folding of the graphene sheet must be carried out in or- 



der to obtain the shell. The shells are usually classified 
into three different types: (n,0) or zig-zag shells, (n,n) 
or arm-chair shells, and chiral shells, of indices (n,m) 
where n > m > 0. Zig-zag and arm-chair shells are said 
to be achiral (they can be superimposed onto their mirror 
images). The ordering of the shells in a MWCNT is usu- 
ally turbostratic, i.e. the pattern of atomic arrangement 
may vary (and in general does vary) from one shell to 
the next, or in other words, different shells usually have 
different chiralitieso. 

After the discovery of MWCNT's, a procedure for sva- 
thesizing single- wall nanotubes (SWCNT's) was foundEl 
These tubes are found to aggregate into bundles or ropes, 
and their diameter distribution peaks at around 1.4 nm, 
although more recently this distribution has-been found 
to vary according to the synthesis conditions^. The pro- 
duction-of SWNT's has allowed the experimental corrob- 
oratior£9 of aptheoretical prediction made by Hamada 
and coworkers^ soon after the discovery of MWCNT's, 
that the electrical conductivity properties of SWCNT's 
are dependent on the (n,m) indices. This prediction 
stated that SWCNT's can be either metallic or semi- 
conducting; if the indices (n,m) of a nanotube obey the 
relation n — m = 3q (q = 0, 1, 2, 3, . . . ), then the tube 
is metallic, otherwise the tube is semi-conducting. This 
dependence of the electric characteristics of SWCNT's 
upon their structure has raised an interest in the possi- 
bility of devising new synthetic methods that allowed a 
structural selection of nanotubes, not only according to 
their diameter, but also to their chirality. Aiirst step in 
this_direction was achieved hxRedlich et aln5, Carroll et 
al.ni, and by Terrones et al.cH, who, by adding a certain 
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amount of boron during the synthesis, obtained boron 
doped MWCNT's which, interestingly, have an improved 
crystallinity and a larger aspect ratio (quotient of length 
to diameter) with respect to tubes obtained in the ab- 
sence of boron. It was shown that boron appears mostly 
in the form of clusters associated with the tips of the nan- 
otubes. But most importantly, a recent combined expex=. 
imental and theoretical study by Blase and coworkersEj 
has demonstrated that the MWCNT's thus obtained con- 
sist mostly of zig-zag shells. This study also sheds some 
light on the role played by boron in favoring the zig-zag 
structure over others, and tries to explain the larger as- 
pect ratio observed in boron assisted synthesis. This lat- 
ter issue, though, was addressed by First-Principles (FP) 
Density Functional Theory (DFT) Molecular Dynamics 
simulations, and the large computational costs of this 
technique prevented Blase and coworkers from pursuing 
a detailed enough study which clarified completely this 
question. 

In this paper we address the problem of boron-assisted 
nanotube growth using a Tight Binding modeO Tight 
Binding (TB) is an approximate method which never- 
theless is capable of providing extremely accurate results 
in favorable systems. Its main advantage with respect to 
FP methodologies is its comparatively low computational 
cost, which often allows a more extensive study than is 
practical or even possible with higher levels of theory. 
In this work we have used the Density-Functioaal Tight 
Binding (DFTB) model due to Porezag et alEJ, about 
which we give more details in Section |n]. This model has 
proved to be very accurate for carbon based systems. We 
have used DFTB to perform a series of static and dynam- 
ical simulations of SWCNT's, with and without boron 
present, with the aim of understanding the effects of the 
presence of boron on the structural properties of the re- 
sulting NT's. The structure of this paper is as follows: 
in Section || we describe briefly the DFTB model, and 
provide previous examples of its successes in order to jus- 
tify its use here. We also describe the calculations which 
are reported in the remaining of the paper. Section III is 
devoted to a discussion of our simulation results, and we 



summarize our conclusions in Section IV 



II. COMPUTATIONAL DETAILS 
A. Model 

DFTB is is a non-orthogonal Tight Binding scheme 
in which a parametrisation is constructed directly from 
DFT calculations using atomic-like orbitals in the ba- 
sis set, and adopting a two-center approximation for 
the Hamiltonian matrix elements. For more details on 
the para metri sation used here the reader should consult 
referencesE^lBil. The DFTB scheme has proved to be ex- 
tremely successful in the modeling of carbon-based sys- 
tems, in particular carbon clusters and nanostructures. 



Fowler et alSn have used it to analyze the energetic or- 
dering of all 426 cage structures contaming 5, 6 and 
7-membered rings in C40. Ayuela et alS3 have found, 
using DFTB and other five semi-empirical methods, a 
heptagon-containing isomer of C62 which was predicted 
to be more stable than any of the other 2385 classi- 
cal fullerene isomers {i.e. isomers containing only pen- 
tagons and hexagons). DFTB has also been used to de- 
termine the mechanical properties of single-wall C, BN 
and some B^CyNz nanotubesEU, providing results which 
are in excellent agreement with the available experimen- 
tal data. More recently DFTB has also been used to 
study the structural, mechanical and electronic prop- 
erties of a novel family of laminar carbon structures 
known as Haeckelites , as well as those of their tubu- 
lar counterpartsO. Haeckelites consist of pentagons and 
heptagons in equal number, with an arbitrary number 
of hexagons. The suitability of DFTB for performing 
Molecular Dynamics simulations in carbon-based sys- 
tems-has been most recently demonstrated by Fugaciu 
et alE3, who have used it to study the conversion of di- 
amond nanoparticles to concentric shell fullerenes. The 
many examples of the use of DFTB in the context of 
carbon based systems and the accuracy of the results 
reported give us confidence in the reliability of this the- 
oretical model. 1 — I 

Like in other Tight Binding modelsca, in DFTB the to- 
tal energy is calculated as the sum of two contributions: 
the band structure contribution, which is mostly attrac- 
tive, and the repulsive pair-potential contribution, which 
accounts for the core-core repulsion and the double- 
counting of the electron-electron interaction which is im- 
plicit in the band structure term. The band structure en- 
ergy is calculated by straight forward diagonalisation of 
the Hamiltonian, summing the eigen-values of the occu- 
pied states weighted according to their occupation num- 
bers, i.e 



(1) 



where /„ is the population of state n (which in our case 
is equal to 1 for occupied states, and to zero for unoccu- 
pied, i.e. we have assumed K electronic temperature), 
and e n is the eigen- value of state n. The factor of two 
accounts for the degeneracy of spin. The band structure 
contribution to the atomic force is then calculated 
from the Hcllmann-Fcynman theorem: 



2f n Ci (ViH-c„ViS)C, 



(2) 



where C„ is the vector representation of eigen-state n, 
H is the Hamiltonian matrix, and S is the overlap. The 
repulsive pair potential energy and force contributions 
are trivially added to equations [I] and |^, respectively. 
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B. Calculations 



A. Carbon nanotubes in the absence of boron 



We have performed two types of calculations: a) struc- 
tural relaxation calculations, in .which, using the Conju- 
gate Gradients (CG) techniqueE£l, the positions of the 
atoms in a system are displaced until a minimum in the 
potential energy hyper-surface is found; and b) molecu- 
lar dynamics (MD) simulations, which we have employed 
to study the time evolution of the systems considered 
here at a range of temperatures. More specifically we 
have performed canonical ensemble molecular dynamics 
calculations, in which the conserved quantities are the 
number of atoms, the volume and the average temper- 
ature (NVT-MD). We have implemented th.e_NVT-MD 
algorithm of Nose as modified by HooverEZfE3. In this 
algorithm the physical system of interest is extended by 
the addition of an extra degree of freedom, playing the 
role of a thermostat, which interacts with the physical 
system in such a way as to fix the average temperature. 
In this algorithm the mass associated with the thermo- 
stat is somewhat arbitrary, as it does not affect the value 
of the average temperature, only the size of the instan- 
taneous temperature fluctuations. We have chosen this 
mass to be equal to the total mass of the physical system. 
We have used the implementation of theJNose-Hoover 
algorithm described by Frenkel and SmitEj. Although 
the total energy is not conserved in NVT-MD, there is a 
conserved quantity which can be monitored to ensure the 
correctness of the implementation. This magnitude plays 
the role of total energy of the extended system. In our 
simulations we have used a time-step of 1 fs, which has 
proved to be sufficiently small to maintain the conserved 
quantity oscillations smaller than 1 part in 10000 during 
the length of the simulations, with no appreciable drift. 

Our studies have covered a range of temperatures, 
namely 1000, 2000, 2500 and 3000 K. At each tempera- 
ture the dynamics of the system were monitored during 
at least 10 ps. The simulations were carried out sequen- 
tially, using the coordinates and velocities at the end of 
a run at a given temperature as a seed for the simulation 
at the next higher temperature. The transition from one 
temperature to another was made by gradually increas- 
ing the temperature of the thermostat at a rate of 0.5 
K/fs. To use a lower heating rate would be prohibitive, in 
view of the computational costs involved. Nevertheless, 
once the thermostat reached the desired temperature, the 
system was allowed to reach equilibrium at the new tem- 
perature during 1 ps, in order to minimize as much as 
possible the effects of such a high heating rate, before 
performing the simulation at the new temperature. 

III. SIMULATION RESULTS 



Before addressing the effect of boron on the structure 
and growth mechanism of nanotubes, it is necessary to 
study the dynamics of carbon nanotubes in the absence of 
boron. In so doing we can also assess the quality and ap- 
propriateness of the model. In order to study the dynam- 
ics of the open ended tubes, we first considered the sys- 
tems illustrated in Figure [j]. These systems are a (10,0) 
nanotube, with one end saturated by ten H atoms, and 
a (5,5) nanotube, also with one end saturated by ten H 
atoms. Both tubes consist of 120 C atoms, plus the 10 
H atoms. Figure [l] shows the structures obtained after 
a CG relaxation, which were then used as starting con- 
figurations for the NVT-MD simulations. To facilitate 
the analysis of the MD simulations, the 10 H atoms were 
kept fixed throughout the dynamical simulations. While 
this is an approximation, we have chosen a tube section 
as large as was practical in order to minimise the possible 
effects of having the H atoms frozen on the dynamics at 
the other end of the tube. 




FIG. 1. The relaxed structures of a (10,0) nanotube (a) 
and a (5,5) nanotube (b), used as starting configurations for 
the Molecular Dynamics simulations described in the text. 

We first discuss the case of the (10,0) nanotube. At the 
lowest temperature considered (1000 K) the structure of 
the nanotube remains unaltered except for the normal 
thermal oscillations around the equilibrium atomic posi- 
tions. All the hexagonal rings present in the initial struc- 
ture (Figure Q) retain their identity, as can be seen in 
Figure 0(a), which illustrates the configuration attained 
after 10 ps of MD run at this temperature. In view of the 
fact that the simulation of the open nanotube end at 1000 
K failed to produce any structural reordering within this 
time scale, we proceeded to heat the system up to 2000 K. 
At this temperature, structural changes begin to appear 
in the open edges of the nanotube, as can be seen in Fig- 
ure ^(b). Indeed, some hexagonal rings at the edge have 
been broken, resulting in chains of carbon atoms, while 
others have fused into pentagon-heptagon (5/7) pairs. 
After 10 ps at 2000 K, two pentagonal rings can be seen, 
as well as one heptagon. These non-hexagonal rings are 
produced by the fusing of two hexagons, to give a pen- 
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tagon and a heptagon. However, we find that heptagons 
appear to be less stable, and break more easily into car- 
bon chains that remain attached to the nanotube edge, 
while the pentagons remain stable. The presence of these 
pentagons induces curvature in the structure, so the edge 
of the nanotube bends inwards. 



(a) (b) 




FIG. 2. (a) Final structure of the (10,0) nanotube obtained 
after 10 ps of dynamics at 1000 K. (b) Final structure of the 
(10,0) nanotube obtained after 10 ps of dynamics at 2000 K. 

Clearly, a tendency to form a fullerene-like cap is ob- 
served, a situation that is after all energetically more 
stable than the open edge. Therefore, we continued the 
simulation at a temperature of 2500 K. Like before, the 
temperature of the thermostat was increased at a rate 
of 0.5 K/fs, and once the desired thermostat tempera- 
ture was achieved the combined system of nanotube and 
thermostat was allowed to equilibrate during 1 ps. The 
dynamics of the system were then monitored during a 
further 20 ps. Figure || illustrates six representative con- 
figurations of the system at this temperature. As can be 
seen, the process of tube closure initiated at 2000 K is 
further facilitated at 2500 K. The chains of carbon atoms 
which were seen already at 2000 K can now form links 
bridging across the open edge. This, combined with the 
presence of the pentagonal rings at or close to the edge, 
has the overall effect of bending inwards the nanotube 
walls. As the simulation proceeds, the number of freely 
oscillating chains is reduced, due to their tendency to 
bond into the dome-like structure being formed. Eventu- 
ally the tube closure is completed, and all carbon atoms 
are three-coordinated in an sp 2 hybridisation fashion. Al- 
though the resulting closed structure is highly strained, 
due to the presence of some small rings (a tetragon can be 
seen) and to the adjacency of pentagons, this is certainly 
a deep local minimum of the potential energy hyper- 
surface, as the structure, once closed, remains unaltered 
(except for thermal oscillations), even after heating up 
the system to 3000 K and following the dynamics at this 
temperature for a further 10 ps. Obviously, the time-scale 
for structural re-orderings once the nanotube has closed 
is much larger than the time-scale for the closure itself. 
This is because structural changes now must happen via 
the Stone- WalesEj mechanism, which has a high activa- 
tion barrier, and therefore occurs very infrequentlytl 



(a) (b) 




(c) (d) 
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(e) (f) 




FIG. 3. Different stages of the closure of the (10,0) nan- 
otube at 2500 K. These structures were observed at approx- 
imately (a) 2 ps, (b) 5 ps, (c) 8 ps, (d) 11 ps, (e) 15 ps and 
(f) 18 ps. 

We now discuss the results obtained from the simu- 
lation of the (5,5) nanotube. During the first 10 ps of 
dynamics, in which the average temperature was fixed 
at 1000 K, one of the edges of a hexagon at the end of 
the nanotube is broken, and as a result a carbon chain is 
formed which oscillates under the influence of the thermal 
motion, but no (5/7) complexes are observed at this tem- 
perature. When the temperature is increased to 2000 K 
two more hexagons break, but still no (5/7) pairs are 
formed. At 2500 K nearly all edge hexagons have broken, 
and consequently there is an abundance of carbon chains. 
Among these the first two pentagons can be observed, but 
there are still no signs of heptagons. This is probably be- 
cause the mechanism of pentagon formation in this case 
is somewhat different than in the (10,0) nanotube. In 
the latter, pentagons and heptagons occur in pairs, re- 
sulting from the fusion of two edge hexagons. Here the 
pentagons appear from the re-bonding of chains result- 
ing from the breaking of edge hexagons. The presence 
of these pentagons induces curvature in the structure, 
which will eventually cause the tube closure. At 3000 K 
a total of three pentagons are seen, and the structure is 
rapidly evolving towards complete closure, but this only 
happens beyond 10 ps of dynamics at this temperature. 

Our results show unambiguously that the open edges 
of nanotubes are unstable, and at sufficiently high tem- 
peratures close spontaneously by forming half-fullcrcnc 
caps. The time scales we find for the nanotube closure 
are only orientative, as they may vary slightly for dif- 
ferent initial conditions, and certainly they are sensitive 
to the temperature, but we can say that above 2000 K, 



closure occurs within a few tens of ps. This finding is 
in agreement with previous first-principles resulta£3. At 
temperatures between 1000 and 2000 K the time scale 
for closure (which remains beyond the scope of the sim- 
ulation times covered here) could be in the order of hun- 
dreds of ps or even in the ns range, but even so this 
seems to us to be too short a time scale for permitting 
nanotube growth, and this is presumably why single-wall 
nanotubes cannot be obtained without using transition 
metal nano-particles which act as catalysts. However, 
the microscopic details of this catalytic growth process 
are as yet far from being understood^ 

In order to gain some insight into the experimental ob- 
servation that boron doping can result in the production 
of longer and more crystalline multi-wall nanotubes, and 
in particular how it favors the zig-zag type tubes over the 
arm-chair or chiral ones, we have performed a number of 
static and dynamical simulations in the (10,0) nanotube 
system with different degrees of boron substitution. 

B. Boron doped nanotubes 

We first address the question as to why boron tends 
to accumulate at the ends, of the nanotubes, as has been 
experimentally observedcS. To investigate this site pref- 
erence, we have performed a series of relaxation calcula- 
tions in which we compare the energy of nanotube and 
fiat edge structures with a boron atom substituting a car- 
bon atom at different sites. We have considered two dif- 
ferent nanotube structures, one a zig-zag nanotube with 
indices (10,0) (the same as illustrated in Figure [j]), and a 
(6,6) arm-chair nanotube. For each of the two structures 
we have taken into account five possible boron substi- 
tutional sites, namely: the boron atom substitutes at an 
edge carbon site (site a), the boron is located at a site one 
bond away from the edge site (site /3), two bonds away 
(site 7), three (site 5) and six bonds away from the edge 
(site e). These positions are illustrated schematically in 
Figure |I| 
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FIG. 4. (a) Schematic view of the boron substitutional sites 
considered in the (10,0) nanotube and the zig-zag graphene 
edge; (b) the substitutional sites in the (6,6) nanotube and 
arm-chair graphene edge. 

Our results, given in Table |, show that in both types of 
nanotubes the most favorable position for boron to sub- 
stitute at is in the vicinity of the edge . However, there 
are marked differences between both types of edges. In 
the zig-zag case, the energy of the structure with boron 
substituted in site (3 is 1.4 eV higher than the structure 
with boron substituted at the a site, almost the same 
as for substitution of boron far (site e) from the edge 
(1.41 eV). So, in the case of the zig-zag edge it is much 
more favorable for boron to substitute at the edge than 
elsewhere in the nanotube. In the arm-chair edge, how- 
ever, it turns out that the most stable structure is that in 
which the boron atom is placed at a site one bond away 
from the edge, i.e. site f3 in Figure |J(b), which turns out 
to be 0.58 eV more stable than the structure obtained by 
placing the boron atom directly at the edge (site a) . Nev- 
ertheless, boron will still preferentially locate itself close 
to the edge rather than far away from it, as the energy 
difference between the most stable configuration and that 
with boron placed in site e is 0.88 eV. These findings can 
be easily rationalized: in the zig-zag case it is more fa- 
vorable for boron to be at the edge than elsewhere in the 
tube because in this way the number of carbon dangling 
bonds at the nanotube edge is reduced, given that boron 
has one electron less than carbon. In the arm-chair case, 
however, the same argument does not apply, because the 
carbon atoms at the edge are able to pair up forming a 
stable double bond. By placing a boron atom at the edge 



a double bond cannot form, and the configuration that 
results is thus less stable than configuration (3. The same 
trends are evident from the calculations involving the flat 
graphene edges, indicating that the curvature of the nan- 
otubes does not play any significant part in determining 
the energetics of the boron-substituted structures. 

TABLE I. Energies of boron-substituted structures for 
zig-zag and arm-chair nanotube and graphene edges. The 
energies are given in eV relative to the most stable structure 
in each case. 



bonds 


(10,0) 


zig-zag- 


(6,6) 


arm-chair 


from edge 


nanotube 


edge 


nanotube 


edge 





0.00 


0.00 


0.58 


0.48 


1 


1.39 


0.99 


0.00 


0.00 


2 


0.57 


0.67 


0.60 


0.66 


3 


1.30 


1.42 


0.83 


0.75 


6 


1.41 


1.35 


0.88 


0.88 
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These static calculations illustrate how boron could 
play a role in favoring zig-zag structures over arm-chair 
ones, as has been found experimentally. Our findings 
are in very good agreement with the DFT plane-wave 
pseudo-potential results reported by Blase and cowork- 
ersc3. However, these results do not give any explanation 
for the increased length of the nanotubes synthesized in 
the presence of boron. Thus, in order to gain some in- 
sight into how boron can assist the production of longer 
nanotubes, we have performed a series of MD simula- 
tions of a (10,0) nanotube with varying degrees of boron 
substitution at the edge. We have employed the same 
structure as illustrated in Figure |l|, but placing four, six, 
eight and ten boron atoms at the tube edge, substitut- 
ing as many carbon atoms in each case. Starting from 
such configurations, we have performed MD simulations 
following the same pattern described above for the un- 
doped carbon nanotubes. Let us comment briefly on the 
details of the dynamics in each case. In all instances, at 
the lowest temperature considered (1000 K), the struc- 
tures remain unaltered during the time spanned at this 
temperature (10 ps); all hexagonal rings present initially 
maintain their identity during this time, and it is not 
until the system is run at a temperature of 2000 K that 
structural changes begin to appear. 

When only four boron atoms are present, after 10 ps 
of dynamics at 2000 K two (5/7) pairs are formed at the 
tube edge. The structure of these (5/7) pairs is such 
that the common edge between both rings always incor- 
porates a boron atom [see Figure ||(a)]. This is a struc- 
tural feature that repeats itself frequently in the other 
cases considered, as will be seen below. At this temper- 
ature, a transient C4 ring, and a B2C6 octagonal ring 
are also seen. When the temperature is raised further to 
2500 K, the boron containing pentagons still survive, but 
the heptagonal rings break, and carbon chains linked to 
the pentagons are formed. These chains, of length equal 
to two or three bonds, oscillate widely, eventually linking 
with similar chains or unsaturated bonds along the open 
edge, thus initiating the tube closure. The tube is closed 
after approximately 8.5 ps at this temperature. Given 
the low concentration of boron at the tip, the dynamics 
observed here is rather similar to that already described 
for the pure carbon case. 



(a) 




FIG. 5. Schematic view of (a) the (5/7) edge complex ob- 
served during the MD runs of the (10,0) nanotube (the circle 
indicates the position where boron is found in the boron doped 
cases); (b) one of the the possible edge reconstructions involv- 
ing pentagons and heptagons in the arm-chair type tubes, the 
(7/5/7) complex, and (c) the (5/7) arm-chair edge complex. 
The circles indicate the substitutional sites considered in the 
relaxation calculations (see text). 

When the number of boron atoms is increased to 6, 
at the temperature of 2000 K two (5/7) complexes like 
those observed in the four-boron atom case are formed, 
but the structure is otherwise unaltered. These boron- 
containing (5/7) pairs are remarkably stable; they remain 
intact even after 10 ps at 2500 K. In this particular case, 
it is necessary to rise the temperature further, to 3000 K, 
for the tube closure to occur. This happens only after 
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8 ps of dynamics at this temperature, which is indicative 
of a higher stability of the open edge when compared to 
the tube with only four boron atoms. In the case of the 
nanotube doped with eight boron atoms, the structure 
develops two (5/7) pairs within 10 ps at 2000 K, much like 
in the previous cases. However, once the temperature is 
elevated to 2500 K, the structure closes completely within 
7.5 ps. This result appears to contradict the thesis that 
increasing the amount of boron at the tube edge delays 
the closure process, but in fact this is not necessarily so, 
as will be argued below. 

Finally, let us consider the case of the tube with ten 
boron atoms. Like in all the cases analyzed earlier, 
structural changes at the open edge begin to appear at 
a temperature of 2000 K, at which a boron-containing 
(5/7) pair is formed. Three such complexes can be seen 
after 10 ps at 2500 K, but the structure is otherwise un- 
changed. It is necessary to reach a temperature of 3000 K 
and monitor the dynamics of the system for 11 ps at 
this temperature before the tube can be said to he. com- 
pletely closed. We note that Blase and coworkerscfl, who 
have also performed MD simulations of this particular 
system using DFT with a basis set of plane-waves and 
the pseudo-potential approximation up to temperatures 
of 2500 K for 10 ps, failed to observe the closure. The 
results reported here indicate that, at such a tempera- 
ture, it is most likely that the closure would not occur 
below 10 ps, as it required 11 ps at 3000 K to observe 
the closure in the simulations that we have performed. 

(a) (b) 




FIG. 6. The dome structures of the closed boron doped 
nanotubes arrived at during the MD simulations: (a) nan- 
otube with four boron atoms, (b) nanotube with six boron 
atoms, (c) nanotube with eight boron atoms and (d) with ten 
boron atoms. 



Figure |6| illustrates the structures obtained at the end 
of the MD simulations described above. As can be seen, 
regardless of the amount of boron present initially at the 
tube edge, the final structures are invariably closed, but 
as seen in the discussion, there are notable differences 
observed during the dynamics in each case. That the 
nanotubes evolve in such a way as to achieve a closed 
structure is, after all, not so surprising, given that the 
closed structures, with all atoms connected in an sp 2 
network are always more stable than the corresponding 
open structures. For example, the structure illustrated 
in Figure ^(d), once relaxed at K, is 13.7 eV more sta- 
ble than the corresponding open structure (a figure to be 
compared with an energy difference of 17 eV, obtained 
for the pure carbon case using the same two structures) . 
Nevertheless, as can be seen from the discussion above, 
the details of the closure process vary according to the 
amount of boron present at the edge. 

Our MD results seem to indicate that a delay of the 
tube closure occurs when boron is present in the sys- 
tem. However, it is difficult to establish this beyond 
doubt purely on the basis of a small number of MD tra- 
jectories, because at each temperature and boron com- 
position there will be a distribution of possible closure 
times (depending on the initial conditions) for a given 
type of nanotube. To get a clear picture of how the clo- 
sure time distribution changes as the amount of boron is 
increased would require a very large sample of MD simu- 
lations, much larger than would be practical to carry out 
given the computational costs involved. Our simulations 
are clearly insufficient in number to reveal unambigu- 
ously the changing trends in closure times with increasing 
boron content, which could explain the apparently con- 
tradicting result obtained for the tube with eight boron 
atoms, which closes faster and at a lower temperature 
than the tube with only six boron atoms. Nevertheless, 
these simulations are sufficient to clarify the role played 
by boron when present at the open edge of a nanotube, 
and are thus extremely revealing. They indicate that be- 
fore the tube closure can begin to occur, (5/7) complexes 
form at the edge. This is true of both the pure carbon 
and the boron doped tubes. Under the effect of the rapid 
thermal motion the heptagons seem to be more prone to 
breaking, and the chains of atoms thus formed initiate or 
can participate in the tube closure. Interestingly, when 
boron is present, it seems to be always associated with 
these structures once they are formed, in the manner de- 
picted in Figure ||(a). 

These observations have lead us to perform a new set 
of static structural relaxation calculations, in which we 
investigate the stability of these edge (5/7) complexes in 
both (n,0) and (n,n) nanotubes. We have compared the 
energy of a normal zig-zag edge structure, such as that 
shown in Figure [j] with that of a zig-zag edge containing 
a (5/7) pair, both in the pure carbon case and in the 
case in which the (5/7) pair contains a boron atom sit- 
uated between the pentagon and heptagon at the edge 
[Figure ||(a)]. It turns out that a (5/7) pair placed at 
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the edge of a pure (10,0) tube lowers the total energy by 
about 0.56 eV. When boron is present between the pen- 
tagon and heptagon, the energy gain is slightly larger, 
0.67 eV, according to our calculations. This seems to in- 
dicate that the formation of such boron containing edge 
complexes stabilizes the open edge, and consequently de- 
lays its closure. 

A (5/7) complex is also possible at an arm-chair nan- 
otube edge, where it is also possible to construct a more 
complicated edge structure, involving a pentagon and two 
heptagons, which we call a (7/5/7) edge complex; both 
structures are illustrated schematically in Figure ||(b) 
and |s|(c) . Structure ||(b) is reminiscp»t of the ring-defect 
structure discussed by Crespi et alE3. We have investi- 
gated the stability of these, both in the pure carbon case 
and when boron is present in the different possible sub- 
stitutional sites. The results from these calculations are 
listed in Table H. It turns out that the only reconstructed 
arm-chair edge which is more stable than the structure /3 
shown in Figure [|(b) is that illustrated in Figure |s|(b), 
with the boron atom placed between the two heptagonal 
rings, but not forming part of the pentagon (i.e. site /?). 
Even so, the energy difference is only small (0.33 eV), 
certainly smaller than the energy gain obtained by the 
reconstruction in the zig-zag edge. 

TABLE II. Energies of boron-substituted reconstructed 
edges for the (6,6) nanotube. The labeling of the edges corre- 
sponds to that of Figure JE] (b) and (c). The energies are given 
in eV relative to the most stable unreconstructed substituted 
arm-chair edge structure. 



Position 


Structure (b) 


Structure (c) 


a 


1.10 


3.40 





-0.33 


3.17 


7 


0.46 


2.67 


QJ 




2.64 



The results obtained from these structural relaxation 
calculations provide a possible mechanism for explaining 
both the observed preference for the zig-zag structure and 
the improved aspect ratio of nanotubes synthesized in the 
presence of boron. The reconstruction of a zig-zag edge 
in which boron is present from the saw-tooth structure 
illustrated in Figure [l] to one containing a (5/7) pair is 
exothermic, resulting in a significant stabilization of the 
nanotube. In the case of an undoped nanotube edge, this 
reconstruction also lowers the energy, but by a smaller 
amount. Note also that, while similar reconstructions 
are possible in the case of a boron doped arm-chair nan- 
otube, only one such reconstructions lowers the energy 
with respect to the unreconstructed edge, and only by a 
small amount. Boron, thus, is seen to significantly sta- 
bilize zig-zag edges compared to arm-chair ones and we 
expect the former tubes to grow better; furthermore the 
edge (5/7) relaxation mechanism we have found is only 
effective in zig-zag edges, a fact that contributes to favor 
this type of tube even further. 

IV. CONCLUSIONS 

We have performed an extensive theoretical study of 
the effects of the presence of varying amounts of boron in 
the edge of growing nanotubes using TB MD and struc- 
tural relaxation calculations. In spite of the approximate 
nature of the TB model used here, it gives results that 
compare quite well-with the available FP data from the 
work of Blase et alE3. 

We have shown that both carbon and boron-doped 
nanotubes spontaneously close in a time scale of a few 
(< 20) ps in a temperature range of 2500-3000 K. In the 
case of boron-doped tubes the MD simulations reveal a 
pattern of longer closure times as the amount of boron 
at the tip of the nanotube is increased. The presence of 
boron in the proximity of the tube edge has a stabilizing 
effect, but this effect is larger in the case of zig-zag nan- 
otubes than in the case of arm-chair ones. Furthermore, 
we have shown that in (n,0) edges a reconstruction of 
the edge involving (5/7) pairs stabilizes further the open 
edge. Although similar reconstructions are topologically 
possible in the case of (n,n) edges they do not result in 
any significant stabilization of this type of edge. 

These results contribute t o , a. , bett er understanding of 
the experimental observationsBEDea that boron accumu- 
lates at the nanotube tips, that it favors zig-zag nan- 
otubes over other structures, and that the nanotubes 
synthesized in the presence of boron have a larger aspect 
ratio than other nanotubes. 

* to whom correspondence should be addressed; emai: 
ehe@icmab.es 
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